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This paper furnishes the basis for a companion one, which discusses the pos- 
sibility of describing material particles as localized oscillatory disturbances in a 
mechanical medium. If a medium is to support such disturbances it must reflect 
a part of the energy of a diverging spherical wave. It is here shown that this 
property is possessed by a medium, such as that proposed by Kelvin, in which 
the elastic forces are of gyrostatic origin. This is due to the fact that, for a 
small constant angular displacement of an element of this medium, the restoring 
torque, instead of being constant, decreases progressively with time. 

Introduction 

IN A companion paper 1 it is pointed out that it may be possible to de- 
scribe the behavior of material particles as that of moving patterns of 
wave motion, provided a medium can be found which is capable of sus- 
taining a localized oscillatory disturbance. In most media this is not possible, 
for the energy of the disturbance would be propagated away in all directions. 
Something special in the way of a medium is therefore called for. It must 
be capable of trapping the wave energy released from a central source. 
Kelvin proposed a mechanical medium, the equations of which, for small 
disturbances, were identical with those of Maxwell for free space. The 
medium derived its elasticity from gyrostats. He recognized that, for finite 
disturbances, the restoring torque depends on the time as well as the angular 
displacement. It is the present purpose to show that this time dependence 
imparts to his medium exactly the energy trapping property required. 

The Gyrostatic Ether 

The concept of an ether with stiffness to rotation originated with Mac- 
Cullagh' 2 in 1839, and was further developed by Kelvin 3 in 1888. MacCullagh 
showed that certain optical phenomena associated with reflection could not 
be represented by the elastic solid ether of Fresnel, but required for their 
mechanical representation a medium in which the potential energy is a func- 
tion of what is now called the curl of the displacement. Fitzgerald 1 remarked 
in 1880 that its equations are identical with those of the electromagnetic 

1 R. V. L. Hartley, Matter, a Mode of Motion — this issue of the Bell System Technical 
Journal. 

2 Collected Works of James MacCullagh, Longmans Green & Co., London, 1880, p. 145. 

3 Mathematical and Phvsical Papers of Sir William Thomson, Vol. Ill, Art. XCIX, 
p. 436, and Art. C. p. 466. 

4 Phil. Trans. 1880, quoted by Larmor, Ether and Matter, Cambridge Univ. Press, 
1900, p. 78. 
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theory of optics developed by Maxwell. This conclusion is confirmed in 
later discussion by Gibbs, 5 Larmor, 4 and Heaviside." 

Kelvin, apparently unaware of MacCullagh's work, was led by similar 
considerations to the same result. He went farther and devised a physical 
model which consisted of a lattice, the points of which were connected by 
extensible, massless, rigid rods in such a manner that the structure as a whole 
was incompressible and non-rigid. Each of these rods supported a pair of 
oppositely rotating gyrostats. By a gyrostat he meant a spinning rotor 
mounted in a gimbal so that it is effectively supported at its center of mass 
and can have its spin axis rotated by a rotation of the mounting. The 
resultant angular momentum of the rotors was the same in all directions. 

This model, considered as a continuous medium, exhibits a stiffness to 
absolute rotation, the nature of which can be described by comparing it 
with the elasticity of a solid. A solid is characterized by a rigidity n such 
that small displacements u, v, w are accompanied by a stress tensor, one 
component of which is 



(dv , du\ 

" U + ay) ■ 



For the ether model the corresponding component is 

(dv du\ „ 
" {to ~ Fy) = 2nie 

where <p is a small angular displacement^ the element about the z axis. 
More generally a small vector rotation Ap is accompanied by a vector re- 
storing torque per unit volume, 

AT = -4«A^. (1) 

The quantity An therefore represents a stiffness to angular displacement 

of the element. 

In the appendix it is shown that the lattice of gyrostats, treated as a 
continuous medium, exhibits this kind of elasticity. It is also shown that 
for infinitesimal displacements, the medium is described by the wave 
equations (8a and 6a). 



**$— S. (2) 

5 Collected Works of J. Willard Gibbs, Longmans Green & Co., New York 1928, Vol. 
' ■ Heaviside, Electromagnetic Theory, Ernest Benn, Ltd., London, 1893, Vol. I, p. 226. 
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where po is the constant density, rjn is a generalized stiffness of the undis- 
turbed medium, given by (7a), q is the vector velocity, and T is the torque 
per unit volume. In a plane wave q is normal to the direction of propagation. 

T 

— is a tractive force per unit area in the direction of q, which acts on a surface 

normal to the direction of propagation. 

If, however, the amplitude is finite the equations become much more 
complicated. For present purposes we need consider only waves for which 
there is no component of velocity or torque in the direction of propagation, 
and we need consider only plane polarized waves for which the direction of 
the velocity is the same at all times and places. Also, as will appear below, 
we are concerned with the equations which describe a wave of infinitesimal 
amplitude which is superposed on a finite disturbance. This description need 
cover only infinitesimal ranges of time and position. It can therefore be 
expressed in terms of wave equations in which the constants of the medium 
have local instantaneous values which depend on the finite disturbance. 

Subject to these restrictions it is shown in the appendix that (2) is to be 
replaced by (23a) 

VX@)-J.p£, (4) 

where l q is a unit vector in the fixed direction of the velocity, and p is an 
instantaneous local density, defined in terms of the finite disturbance by 
(20a). And, in place of (3), (22a) 

where l? is a unit vector in the direction of the axis of rotation, p is again 
an instantaneous local density, c is an instantaneous local velocity derived 
in the usual way from p and an instantaneous local stiffness 77, while /is a 
function defined by the relation, (13a), 

T = -iMr, ')■ 

This function takes account of the fact that when the spin axis of the rotor 
is given a constant finite displacement, the restoring torque is not constant 
as in (1), but changes with time as the spin axis rotates toward the axis of 
displacement, and so reduces the component of the spin which is normal 

df 

to the displacement axis and so is effective in producing stiffness. — 4 — 

dt 

represents the rate of this change in torque for a fixed angular displacement. 

— 4 — is to be interpreted as the rate of change of torque with angular 

d<p 
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displacement, when the time consumed is infinitesimal, that is when the 
angular velocity is infinite. It is therefore an instantaneous local angular 
stiffness from which the instantaneous local generalized stiffness r\ is derived 
as in (19a). 

To simplify these expressions, let the direction of propagation be x and 
that of q be y. Then 

V X q = i — (jq) =^. 
dx dx 

so l t is in the direction of s, and represents a clockwise rotation about z. 
(5) then becomes the scaler equation 

£i = _ J_ 
d.v pc' 2 

T is also in the s direction, so 



te©+'S- 



** G) -«£(»!) --'±G 



But q is in the y direction, so 



dx\2 H dt 



(7) 



These, then, are the desired equations of motion, for the type of wave 
under consideration. 

The Generation of Reflected Waves 

In this section we shall show that when a finite wave is propagated in 
this medium each element of the medium becomes the source of auxiliary 
waves which propagate in both directions from the source. 

To do this we shall make use of the argument by which Riemann 7 showed 
that this does not occur for sound waves in an ideal gas. This will first be 
restated in more modern language. We consider a plane wave propagating 
along the x axis. We picture the finite pressure p and the longitudinal 
velocity u at a point in the medium as having been built up by the successive 
superposition of waves of infinitesimal amplitude, each propagating relative 
to the medium in its condition at the time of its superposition. If the first 
increment is propagating in the positive direction, 

A d P 

au = — , 
pc 

7 Lamb, Hvdrodvnamics, Sixth Edition, p. 481. Rayleigh, Theory of Sound, Second 
Edition, Vol. II. p.'38. 
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where the characteristic resistance is pc. Here 

'"* 

He assumes adiabatic expansion, so that p and c are functions of p only. If 
a second incremental wave of pressure dp, also traveling in the positive 
direction, be added, its velocity increment, being relative to the medium, 
will add to that already present. Its value will be related to dp through a 
new characteristic resistance corresponding to the modified density result- 
ing from the previous increment. Hence the velocity u resulting from a 
large number of such waves will be 



r dp 

= I — = w, 

Jo oc 



Jo pc 

where w is the quantity represented by co in Lamb's version. If, then, all 
of the wave propagation is in the positive direction 

it = IV. 

Similarly, if an incremental wave is traveling in the negative direction, 

* = =*t , 

pc 
and the condition for all the propagation to be in that direction is 

it = —w. 

Obviously, then, if u has some other value than one of these it results from 
the addition of increments some of which propagate in each direction. 
Riemann deduces from the aerodynamic equations that 



Q t + (it + c) £) (to + u) = 0, (8) 

W - c) -f) (w - «) = 0, (9) 



dt dx. 

That is, the value of w + it is propagated in the positive direction with a 
velocity of c -f- it and that of iv — u, in the negative direction with a velocity 
c — u. If, over a finite range of x, a disturbance be set up such that neither 
of these quantities is zero, it must be made up of incremental waves in both 
directions. However, as w + u propagates positively it will be accompanied 
at any instant by a value of w — it which has been propagated from the other 
direction. But, since the value of this was initially finite over a limited dis- 
tance only, when all of this finite range is passed, w — u will be zero, u will 
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be equal to w and all of the wave will be traveling positively. A similar 
argument applies at the negative side of the wave. Thus the initial disturb- 
ance breaks up into two parts which travel in opposite directions without 
reflection. More generally, these considerations hold for any medium in 
which the stress is a function of the strain only. 

For the ether model, since we have assumed the displacements are normal 
to the direction of propagation, the velocity of wave propagation relative to 
the medium is the same as that relative to the axes. 

If now, following Riemann, we let 

(10) 









dw 


pc \ 


I)> 




so that now 






w = 


J pc 


© 


7 


then from (7) 


and 


(6) 


d( i - . 

dt 


dw 
dx ' 












dw 
dt 


dx 


2 
pc 


9f 

dt 



Adding and subtracting gives 



- + c — ) (w + q) = - — rj- , 
dt dx) pc dt 

(d d \ , , 2 df 

\d- t - c ir x ) {w - q) = -7cTt> 



df : . 



which are to be compared with (8) and (9). Hence when -i is not zero the 

values of w + g and w - q are not propagated without change. 

To show that reflection occurs, consider a disturbance at a point x at 
time L characterized by q and w. At x and t + A/, w + q will differ from the 

d 2 df 

value it had at x - cAt, t, or w + q - ~ (w + g)cM, by - — - A/. The 

increment at x in time A/ is 

A ? r) f 

Aw + Ao = - — (w + g)cAl -f- At, 

dx pc at 
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and 

Aw — Ao = — (w — q)cAt — A/. 

dz pc 3/ 

From which 

dx pc dt 

A ^ A# 

Ao = — c — At. 

a / 

Hence the velocity is the same as when -^ is zero but w is changed by 

at 

1 df 

— -At. But the only way in which w can change with q constant is 

pc dt 

by adding waves of equal amplitude propagating in opposite directions, so 

that their contributions to w are equal and those to q are equal and opposite. 

From (10) this involves an increment of — of —2 -^ A/ or a time rate of change 

2 dt 

of —2 — .This agrees with (6), from which it is evident that the presence of 

— alters — from what it would otherwise be by — — - — . But, since q is 
dt dx J pc 2 dt * 

unchanged, the velocities at * + — and x — are increased by — — - — Ax 

2 I pc- at 

1 df 

and — - — Ax. The first is the velocity associated with an auxiliary wave which 
pc 1 dt 

propagates in the positive direction of x, and the second that of one which 

propagates in the negative direction, that is a reflected wave. Hence the 

1 df 

medium generates a reflected wave of — - — per unit length in the direction 

pc 1 dt 

of propagation. 

The Reflection of a Progressive Diverging Wave 

So far attention has been confined to a single point. If a continuous dis- 
turbance is being propagated, it is important to know how the waves reflected 
at different points combine, for it is conceivable that they may interfere 
destructively. From the standpoint of the application to be made of these 
results in a companion paper, the case of most interest is that in which energy 
is propagated outward from a central generator as a sinusoidal wave of 
finite amplitude, beginning at time zero. Near the center, the wave of dis- 
placement will include radial as well as tangential components. As the radius 
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increases the radial components become relatively negligible. We shall 
confine our attention to this outer region, where, in the absence of reflection, 
the propagation differs from that of a plane wave only in that the amplitude 
varies inversely as the radius. We shall neglect the effect of any reflections 
on the outgoing wave, and calculate the resultant reflected wave at a radius 
r\ as a function of the time and so of the radial distance r the wave front 
has traveled. 

If the outgoing wave were of infinitesimal amplitude, its velocity qo 
could be represented by 

qo = — Qo sin (col — kr), (11) 

r 

for values of r < ct, and by zero for r > ct, where Q is the amplitude at 
some reference radius r . The sine function is chosen to avoid the necessity 
of an infinite acceleration at the wave front, as would be required by a 
cosine function. When the amplitude is finite this wave suffers distortion 

due to the fact that k which is equal to - varies slightly with the variations 

c 

in the instantaneous value of c. However, these will be small and, since 
fluctuations in velocity alone do not cause reflection, we shall neglect them. 
The procedure is to make use of q Q to calculate the reflected wave incre- 
ment generated in a length Ar' at a radius r' , calculate the amplitude and 
phase of this at a fixed point r\ <r', and at r\ integrate the waves received 
there for values of / from ri to the farthest point from which reflected waves 
can reach r\ at the time t under consideration. 

To find the reflected wave generated in a length Ar' at r', we have from 
above that its velocity 

pc 2 dt 
From (21a), (19a) and (17a) 

1 Fi 



2 — / r r ~12\ > 

pc ' ' ■ 



ilo 1 - 



a[f V dtT) 



where tj and a are constants of the medium given by (7a) and (15a). From 
(18a) 



df 2 

m —ario<P 
at 



f<pdt, 



so 
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aF\<p~ I <p dt 



1 - 



a\ I <p dt \ 



which reduces to 



d 9' 2 f ,. 

-,= -a<p Jvdt, 



if we neglect second powers of the variables compared with unity. 
To the same accuracy, from (14a) 

From (11) 

dfc m _roQo r cos (aj/ _ kfl) + 1 sin ((tf/ _ ^1 _ 

Here & is 27r over the wavelength so, if as we have assumed r t , and therefore 
also r', is large compared with the wavelength, we may neglect the second 
term. Then 



/ 



*- -W«ln(«l-*rO. 

2cr 

ipdt = ?%cos(«/ - Jfer'), 

2ca>r 



^ = -A (^Y sin 2 (cot - *r') cos (a,/ - *r0, 
</r 8a> \ cr / 

= _A A*Q?Y [cos («/ - kr') + cos 3(co/ - kr% 
8u\cr / 

This, when multiplied by Ar', gives the value at r' of the wave, generated 

in the interval Ar', which propagates in the negative direction of r. This is 

made up of components of frequency to and 3oj. We are primarily interested, 

from the stand-point of reflection, in that of frequency o>, so we shall confine 

our attention to this component, with the understanding that the other 

can be treated in exactly the same fashion. As the fundamental component 

r' 
propagates inward to f\ it increases in amplitude in the ratio — and suffers 

a phase lag of k(r' — n). If we call the resultant of all the reflected waves at 
f*i, q[ , then the contribution to q[ of the wave generated at r' is 
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Aq [ = -± (*&Y * cos («/ + kn - 2kr')Ar'. 
8co \ c ) r x r 2 

This is to be integrated from r x to the farthest point from which a reflected 
wave has reached r± at the instant t under consideration. This point is at 
Jfa + ct). So 

q [ = - «L (5>@?) / 1 cos («/ + ft* - 2kr') dr'. 

8fi CO \ C / Jrj r 2 

Here the integrand is a function of r' and / and the upper limit of integration 
is also a function of /. We therefore make use of the relation 8 

i />• « * = f (£/<*■ ■>) * + ** ■> I - /<«■ «> t ■ 

Putting / for a, r' for x we have 

4 = « (2&Y [p"' I sin M + *, _ 2 *,) _ s i I 

<// 8rx \ c / L^r! r 2 « (fk + c0 2 J 

which, upon integration becomes, 

/ ' / - ° ^ ro0 °V /"i rifa + fen) - 2*[5t(«/ - kn) - Si(2kn)] 



dt 8fi \ c / \fi 

sin (co/ + M) - [Ci(«* + /tfi) - C/(2£r,)] 

2c 



cos (co/ + £ri) — 



w(fi + rf)7 " 



Since c/i is zero when / is — , its value at / will be found by integrating from 

— to /, so 
c 

„; = • (5&Y (_ tos (W _ W + *L _ 2 *r, 

• I co / .S'i(co/ + *rj sin (cot + *fi) <// + Si(2kr i ) 

• [cos (cot + £r0 — cos 2kri] — co I Ci(cot + kfi) cos (co/ + £ri) dt 

+ Ci(2M) [sin(co/ + kr{) - sin 2£rJ J . 



which reduces to 

8 Byerly, Integral Calculus, second edition p. 99. 
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q [ = • (*SY (cos (a,/ - *n) - -*3- + 2k ri 

• [-[S*(«J + kn) - Si(2krJ\ cos (co* + kn) 
- [Ci(w* + kn) - Ci{2kn)] sin (orf + kn) 

+ 5i(2cof + 2kn) - 5«(4*fi)] J . 

The first term represents the value at r x of an outwardly moving wave in 
phase quadrature with the main wave. The second is a transient, the value 
of which is equal and opposite to that of the first term at the instant that 
the main wave passes t\ . The first two terms in the inner bracket are waves 
which propagate inward and so are to be regarded as reflections of the 
main wave. The last two terms represent a velocity which is zero when the 
main wave passes T\ , and subsequently oscillates about and approaches 

- — Si(4&ri). Physically it appears to result from the particular form chosen 

— 

for the main wave, which starts abruptly as a sine wave. The time integral 
of the impressed force, and so the applied momentum, has a component in 
one direction. Presumably if the main wave built up gradually these terms 
would be absent. 

Returning to the reflected waves, their amplitudes are zero when the 
main wave passes n , after which they become finite. Si(x) and Ci(x) os- 
cillate about and approach - and zero respectively as x approaches infinity. 
Hence, as / increases indefinitely, the amplitudes of the reflected waves 
approach - — Si(2kr{) and Ci(2kr{). For the assumed large values of 2kr Y 
these quantities are small compared with unity. When multiplied by 2kri 
their variation is very slow. Hence the amplitudes vary roughly as i , 

and approach zero as the main wave at r x approaches an ideal plane one. 

However, the significant fact is not that the reflected waves are small 
but that they are of finite magnitude. Because of this the main wave will 
not behave exactly as we assumed above, but will decrease slightly more 
rapidly with increasing radius. This should increase the reflection slightly, 
for the existence of the reflected wave is dependent on the decrease in am- 
plitude with distance when the radius of curvature is finite. 

To describe exactly what happens when the generator begins sending out 
waves from a central point would be hopelessly complicated, but we may 
form a general picture. In the early stages where the curvature is consider- 
able, the reflected waves would be quite large and the main wave would be 
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correspondingly attenuated. The arrival of the reflected waves at the gen- 
erator adds a reactive component to the impedance of the medium, as seen 
from the generator, which reduces the power delivered to the medium. 
Meanwhile energy is being stored as standing waves in the medium and 
the rate of flow of energy in the wavefront is decreasing. The energy in 
successive shells of equal radial thickness decreases with increasing r, in- 
stead of being uniform as it would be in the absence of reflection. In the 
limit it approaches zero, but as the rate of decrease depends on the curva- 
ture, the rate of approach also approaches zero. As the rate at which energy 
is stored and that at which it is carried outward at the wavefront both 
approach zero, the resistance which the medium offers to the generator 
approaches zero, and its impedance approaches a pure reactance. 

The total energy stored in the medium depends on how the over-all at- 
tenuation of the main wave is related to its amplitude. If there were no 
attenuation, the impedance would remain a pure resistance, the energy in 
successive shells would all be the same, and the total energy would increase 
linearly with r, and so with the time, and approach infinity. If the attenua- 
tion were independent of r, the total energy would approach a finite value. 
The present case is intermediate between these, the attenuation being finite 
but approaching zero with increasing r. If we assume it to vary as some 
power of the amplitude of the velocity, then W. R. Bennett has shown that 
if this power is less than the first the total energy approaches a finite value. 
If it is equal to the first, the energy approaches infinity as log r, and if it is 
greater than this, the power approaches infinity more rapidly. Until more is 
known as to the actual variation of amplitude with distance, nothing 
definite can be said about the limit of the total energy. 

APPENDIX: EQUATIONS OF THE KELVIN ETHER 

We are concerned with the wave properties of the model for wavelengths 
long enough compared with the lattice constant so that it may be regarded 
as a continuous medium. Its density is equal to the average mass of the 
gyrostats per unit volume. Its elastic properties are to be derived from the 
resultant of the responses of the individual gyrostats. 

We shall therefore begin by considering the behavior of a single element, 
which is shown schematically in Fig. 1. Here the outer ring of the gimbal, 
which is rigidly connected with the lattice, lies in the * y plane. The axis 
about which the inner ring rotates is in the x direction, and the spin axis C 
of the rotor is in the z direction. We wish to examine the effect of a small 
angular displacement <p of the lattice, that is, of the outer ring. If it is about 
x or z, it will, because of the frictionless bearings, make no change in the 
rotor. If it is about y it will produce an equal displacement of the spin axis 
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Fig. 1 — Diagram of a gyrostat, showing its axes of rotation. 

C about y. To study its effect we make use of Euler's equations for a rotating 
rigid body. 9 

A d Jp - (b - C)mm - L, 

at 



B a J^l _ ( C - A)m<* = M, 

at 



C^ - (A -B)wN, 

at 

where a>i , w 2 and oj 3 are the angular velocities about three principal axes of 
inertia, fixed in the rotor, the moments of inertia about which are A, B 
and C, and L, M, and .V are the accompanying torques about the three axes. 
They are also at any instant the values of the torques about that set of 
axes, fixed in space, which, at the instant, coincide with the axes 1, 2, 3, 
which are fixed relative to the body. We let the 3 axis coincide with the 
spin axis C. We choose as the 1 and 2 axes, lines in the rotor which, at the 
instant, are in the x and y directions respectively. Since the moments of 

9 Jeans, Theoretical Mechanics, Ginn and Co.. p. 308. 
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inertia about these are equal, A and B are equal. By virtue of the frictionless 
bearings the external torques L and N about 1 and 3 are zero. 
Introducing these relations we have 

A-£+(C- 4)« 2 W3 = 0, (la) 

at 

A ^T- (C-^)«i«»- M, (2a) 

at 

C £-a (3a) 

From (3a) the velocity of spin o>3 remains constant. The torque M about y 
is then to be found from (la) and (2a). For very small displacements, 

C02 = <p- 

Putting this in (la) and integrating from zero to /, assuming ip to be zero 
at / = 0, gives 

C - A 

0)1 = OJ 3 <p. 



(2a) then becomes 



.... (C - A) 2 2 ,, 

A<p -f- (t) 3 <p = M. 



//-t j \2 2 

This represents an angular inertia A and stiffness — - — -. The system 

will therefore resonate at a frequency — — - . If the frequencies in- 

volved in the variation of <p are small compared with this, the inertia torque 
will be negligible, and the system will behave as a stiffness. If the displace- 
ments about A associated with o>i are very small the restoring torque M 
will act substantially about the y axis. That is, the lattice will encounter a 
stiffness to rotation. 

Since the large number of gyrostats in an element of the model are oriented 
in all directions, an angular displacement of the lattice about y will gen- 
erally not be about the B axis for each gyrostat. If it makes an angle a with 
this axis, then only the component <p cos a of the angular displacement will 
be transmitted to the rotor. The resulting torque will then be S cos a, where 

. _ (C - A) 2 col 
* A • 

It will be directed about B and so will not be parallel to the applied dis- 
placement. However, if a second gyrostat has the position which the first 
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would have if it were rotated about y through w, its torque along y is the 
same as that of the first, and that normal to it is equal and opposite. Hence, 
if the gyrostats are properly oriented, the resultant torque will be parallel 
to the displacement and the medium will be isotropic. The y component of 
the opposing torque will be S<p cos 2 a. Thus if the B axes are uniformly dis- 
tributed in space the total torque will be one third what it would be if they 
were all parallel to the axis of the applied displacement. Hence if there are 
N gyrostats per unit volume the vector restoring torque T per unit volume 
will be 

77 A' (C - /1) 2 W3 _ ,. s 

T- -3 - -j- -<p. (4a) 

The next step is to derive the wave equations for a medium having this 
stiffness to rotation. If the vector velocity q is very small, 

VXj = 2|, (5a) 

where <p is a vector angular displacement of an element of the medium at 
the point under consideration. 2<p plays a role analagous with that of the 
dilatation in compressional waves. Then, from (4a) and (5a), 

where the generalized stiffness of the undisturbed medium, 

N (C - Af 2 
Vo = 12 A C ° 3 " ' ' 

To get the companion equation, we interpret the torque exerted by an 
element in terms of the forces it exerts on the surfaces of neighboring ele- 
ments. Let the x axis Fig. 2 be in the direction of the torque TAx 3 which is 
exerted by the medium within the small cube. This very small torque can 
be resolved into the sum of two couples, one consisting of an upward force 
F y Ax 2 on the right face and an equal downward force on the left one, and 
the other of a leftward force F z Ax 2 on the upper surface and a rightward one 
on the lower one. But, if there is not to be a shearing stress, F v and F z must 

T 
be equal, and each equal to — . Thus a torque per unit volume T is equivalent 

T 
to a set of tangential surface forces per unit area of — each. 

Now consider the force exerted on an element by its neighbors, through 
the adjoining surfaces. To take the simplest case, let T in Fig. 2 be every- 
where in the x direction and independent of z but varying with y. Then 
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the forces exerted on the upper and lower surfaces are equal and opposite. 
That downward on the right face exceeds that upward on the left by 

— ( - Ax 2 ) Ay, so the force in the z direction is — — f - J Ax 3 . By ex- 
tending the argument to three dimensions it is easily shown that the total 



Fy&X' 



F 7 AX 2 -» 




x Fy * x ' 

Fig. 2 — Diagram showing the forces exerted by an element of the medium through 
its surfaces. 



force isV X ( tt 
poAx 3 — , so 



UW.Ifpoi 



is the density of the medium this force must equal 



(?)- 



»x|,|-„J, 



which, since q is small, reduces to 

'X© 



- po dI 



(8a) 



From this and (6a) the velocity of propagation is (vo/po) 1 2 and the char- 
acteristic resistance is (p r] ) 112 . In a plane wave the displacement is normal 
to the direction of propagation. The stress is a tractive force per unit area 

T 

— acting in a surface normal to the direction of propagation. It is in the 

direction of the velocity and in phase with it. 
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However, we are also interested in the case where the amplitudes are 
not negligible. We shall confine our attention to those cases where, as in 
plane or spherical waves at a distance from the source, the velocity is nor- 
mal to the direction of propagation and the variations in the plane of the 
wave front are negligible. (5a) then becomes much more complicated. 

V X q is, however, still a function of — , say 2Fi ( — ] . Then, for small varia- 

at \dt/ 

tions of — in the neighborhood of a particular value, we may write 
at 



-*®2 « 



where F\\ — ) is a function of the particular value of — - . This relation is 
\at / at 

to take the place of (5a). Similarly, if 

**© = *© 

then, in place of (8a), we are to use, for small variations, 

When we come to the transition from (5a) to (6a), however, the situation 
is somewhat different. To see how this comes about, we go back to the 
behavior of the single gyrostat of Fig. 1. It was assumed above that the B 
axis coincided with the y axis However, when the displacement of the 
rotor about A is finite, this is no longer exactly true. The situation is then 
as shown in Fig. 3. A rotation <p of the lattice about y displaces A in the x z 
plane by <p. The accompanying rotation of the rotor about A causes B to 
make an angle with y, which is independent of <p. Then 



From (la) 



Wi = — 



A 

Also 



dip 

OJ2 = -T- C0S V. 

at 



C - A [d^ ... 

a>3 / -T- cos 8 dt. 



6 = I ui dl, 



C - A 



a>3 1 1 -j- cos ddtdt, (11a) 
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which determines 6 as a function of <p and t. From (2a), neglecting the first 
term as above, 

M = S j ^ cos 6 dl, 

and the restoring torque about y, or 

T v = -S cos 6 f ^ cos 6 dt. (12a) 

This, together with (11a), determines T v as a function of <p and /, instead of 
<P alone as it is for infinitesimal displacements. 




Fig. 3 — Diagram showing the displacement of the axes of a gyrostat. 

We assumed here that, in the rest position of the rotor, its B axis coin- 
cides with that of the applied displacement <p. When this is not the case, the 
relations are more complicated, but they should be qualitatively the same. 
Hence, for an element of the medium, the torque per unit volume should 
be a function of tp and / similar to T v , which reduces to — A^ <p for very small 
displacements. Since the restoring torque is in the direction of <p we may 
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write 



T = -IM<P, 



(13a) 



where l v is a unit vector in the direction of the axis of rotation. 

The derivation of the wave equation is much simpler if we consider only 
the case of present interest where the direction of the rotation is everywhere 
the same so that l^ is constant. Then (9a) can be written as 



and (13a) as 



- - '-" -- tfr)ff- 



T = -4/fo I). 



(14a) 



We wish now to replace — by — ( •=■ J • These partial derivatives refer to a 
at at \l / 

constant position so we are interested in the total time derivatives of T as 

given by (12a). To get the desired relation we need to express T explicitly 

in terms of <p and /, that is, we must evaluate <p. Since the variables are 

small, we neglect their products of higher order than the third. Then 



cos = 1 — - a 



f <pdt\, 



where 



Putting 



C - A 
A 



wa 



(15a) 



T = — 4ijo cos 



[ $£ 
J dl 



cos 6 dt. 



in accordance with (12a) and substituting for cos 6 gives 



T = — 4tj \ <p — a<p 



\ <p — atp J ipdt + a I <p 2 I I <pdtj dt \. 



Then 



£--■'-11' 



I (l - a f <pdt j^Z - a<P 2 f tpdt\. 



When ip is constant the first term is zero, so the second term can be inter- 
preted as the partial derivative of T with respect to /. Physically this de- 
scribes the change in torque for a fixed displacement which results from the 
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fact that, as the axis of the rotor rotates toward that of the applied torque, 
the component of the spin which is normal to the axis of displacement pro- 

d<p . 
gressively diminishes. To interpret the hrst term, we let — increase in- 
definitely. The second term then becomes negligible, and when we divide 

through by — , the left side becomes — . But the time increment which 
dt d(p 

accompanies a finite increment of cp is now infinitesimal, and so this may be 

called the partial with respect to <p, with / constant. 

We have then 



dt 
where 



"--»&*+ a 



-£ = —ar}o<p I tpdt. (18a) 

at J 



Substituting for -f from (16a) in (14a), 
at 

d<p 

We may interpret— as an instantaneous stiffness to rotation and define 
dip 

an instantaneous local generalized stiffness by the relation 



dip. 



V = ^7 (19a) 



Similarly from (10a) we may define an instantaneous density by the relation 

p = F' t . (20a) 

Then we may speak of an instantaneous velocity c given by 



c 2 = V - 



(21a) 
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and an instantaneous characteristic resistance pc. Then 

(10a) becomes 

vx (!) = '»>!?> (23a) 

where l q is a unit vector in the fixed direction of the velocity. These are the 
equations of motion which apply to a very small disturbance superposed 
on a finite disturbance. 



